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Abstract 



We propose a new formula for extracting the polarized gluon distribution Ag in the proton 
from the large-pr light hadron pair production in deep inelastic semi-inclusive reactions. Though 
the process dominantly occurs via photon-gluon fusion(PGF) and QCD Compton, we can remove 
an effect of QCD Compton from the combined cross section of light hadron pair production 
by using symmetry relation among fragmentation functions and thus, rather clearly extract 
information of Ag. 



^Talk given by T. Morii at the EMI2001 (International Symposium on Electromagnetic Interactions in Nuclear 
and Hadron Physics), RCNP, Osaka, Dec. 4-7, 2001. 



1 Introduction 



Since the measurement of the polarized structure function of the proton g\ for wide kinematic 
range of Bjorken x and Q 2 by the EMC in 1988 [jTJ , the so-called proton spin puzzle has been still 
one of the most challenging topics in particle and nuclear physics. After the EMC experiment, 
much progress has been attained experimentally and theoretically. [[J A large amount of data on 
polarized structure functions of proton, neutron and deuteron were accumulated. The progress 
in the data precision is also remarkable. Furthermore, with much theoretical effort based on 
the next-to-leading order QCD analyses of the longitudinally polarized structure function gi, 
a number of excellent parameterization models of the polarized parton distribution functions 
(pol-PDFs) were obtained from fitting to many/precise data with various targets for polarized 
deep inelastic scattering (pol-DIS) 0, J|, j[| All of these parameterization models tell us that 
quarks carry a little of the proton spin, i.e. only 30% or smaller. Due to an intense experimental 
and theoretical effort, a rather good knowledge of the polarized valence quark distribution has 
been obtained so far. However, knowledge of Ag is still very limited. As is well-known, the 
next-to-leading order QCD analyses on gi bring about information on the first moment of the 
polarized gluon Ag.jjTl] However, there are large uncertainties in Ag extracted from g± alone. To 
solve the so-called spin puzzle, it is very important to precisely know how the gluon polarizes in 
the nucleon. 

So far, a number of interesting processes such as direct prompt photon production in polarized 
proton-polarized proton collisions open charm|J and J/^HOj productions in polarized lepton 
scattering off polarized nucleon targets, were proposed for studying longitudinally polarized gluon 
distributions. Recently, HERMES group at DESY reported the first measurement of the polarized 
gluon distribution from di jet analysis of semi-inclusive processes in pol-DIS, though only one 



data point was given as a function of Bjorken x. [11] In general, a large-p^ hadron pair is produced 
via photon-gluon fusion (PGF) and QCD Compton at the lowest order of QCD. The PGF gives us 
a direct information on the polarized gluon distribution in the nucleon, whereas QCD Compton 
is background to the signal process for extracting Ag.O, O] For the case of heavy hadron 



pair(DZ), D*D*, etc.) production, we could safely neglect the contribution of QCD Compton 
because the heavy quark content in the proton is extremely small and also the fragmenting 
probability of light quarks to heavy hadrons is very small. However, in this case the cross section 
itself is small at the energy in the running(HERMES) or forthcoming(COMPASS) experiments 



and thus, we could not have enough data for carrying out detailed analysis. For the case of 
light hadron pair production, we will have a rather large number of events because of large cross 
sections. However, in this case QCD Compton can no longer be neglected and hence, it looks 
rather difficult to unambiguously extract the behavior of Ag from those processes. 

In this work, to overcome these difficulties we propose a new formula for clearly extracting 
the polarized gluon distribution from the light hadron pair production of pol-DIS by removing 
the QCD Compton component from the cross section. As is well-known, the cross section of the 
hadron pair production being semi-inclusive process, can be calculated based on the parton model 
with various fragmentation functions. Then, by using symmetry relations among fragmentation 
functions and taking an appropriate combination of various hadron pair production processes, 
we can remove the contribution of QCD Compton from the cross section and thus, get clear 
information of Ag from those processes. Here, to show how to do this practically, we consider 
the light hadron pair production with large transverse momentum. 

2 A new formula for extracting Ag(x) from large— pr pion 
pair production 

Let us consider the process of E+N — > £' + hi + h,2 + X in polarized lepton scattering off polarized 
nucleon targets, where h\ and hi denote light hadrons in a pair. As mentioned above, the spin- 
dependent differential cross section at the leading order of QCD can be given by the sum of the 
PGF process and QCD Compton as follows, 

dAa = dAa PGF + dAa QC D , (1) 

with 

dAa = da-+ — da ++ + <icr + _ — da ++ , 

where, for example, da (_ denote that the helicity of an initial lepton and the one of a target 

proton is negative and positive, respectively. Each term on the right hand side of eq.(P is given 
by 

dAapQp 

~ Ag( V ,Q 2 )dAa PGF £ e^{D^(z 1: Q 2 )D^(z 2 , Q 2 ) + (1 2)}, (2) 

i=u,d,s,u,d,s 

dAa QC D 



~ £ e 2 Af q (r ] ,Q 2 )dAa QCD {D h q ^z 1 ,Q 2 )D h g \z 2 ,Q 2 ) + {l^2)}, (3) 

q=u,d,s,u,d,s 

where Ag(r), Q 2 ), Af q (r), Q 2 ) and D\{z, Q 2 ) denote the polarized gluon and the quark distribution 
function of q with momentum fraction i] and the fragmentation function of a hadron hi with 
momentum fraction z% emitted from a parton i, respectively. dA&pcF and dAcfqcD are the 
polarized differential cross sections of hard scattering subprocesses for ig —>■ i'qq and I ( q — > £' 9 q 
at the leading order QCD, respectively. 

Here we consider the following 4 pairs of the produced hadrons hi with z\ and hi with z 2 , 

(i) (7r + ,7r") , (ii) (7r~,7r + ) , (iii) (it + ,it + ) , (iv) , 

where (particle 1, particle 2) corresponds to (hi with zx, h 2 with z 2 ). Then, the differential cross 
section of eq.(^) for each pair can be written as 

(i) 

dAa T+n ~ 

~ Ag( V , Q 2 )dAa PGF { ^Dl + ( Zl , Q 2 )Df (z 2 , Q 2 ) + Id? ( Zl , Q 2 )Df (z 2 , Q 2 ) 
+ ^Df(z 1 ,Q 2 )Df(z 2 ,Q 2 ) + (7r + (z 1 )^n-(z 2 )) } 

+^Au(rj, Q 2 )dAa QCD {Df ( Zl , Q 2 )Df (z 2 , Q 2 ) + Df (z 2 , Q 2 )Df ( Zl , Q 2 ) } 

+ (contributions from Ad, As, Au, Ad and As) , (4) 

(ii) 

dAa T ' n+ 

~ Ag( V , Q 2 )dAa PGF { ^ {z x , Q 2 )Df(z 2 , Q 2 ) + ±Df {z x , Q 2 )D} + (z 2 , Q 2 ) 
+ lDf(z 1 ,Q 2 )Df(z 2 ,Q 2 ) + (n-(z 1 )^n + (z 2 )) } 

+^Au(rj, Q 2 )dAa QCD [Df [z x , Q 2 )Df(z 2 , Q 2 ) + (z 2 , Q 2 )Df ( Zl , Q 2 ) } 
+ (contributions from Ad, As, Am, Ac? and As) , (5) 
(iii) 

dAa n+7r+ 

~ Ag(r j ,Q 2 )dAa PGF {^DZ + (zx,Q 2 )Df(z 2 ,Q 2 ) + ±Df(z lt Q 2 )Df(z 2 , Q 2 ) 
+ l -D: + { Zl ,Q 2 )Df{z 2 ,Q 2 ) + (n + (zx) <-> n + {z 2 )) } 

+^Au( V , Q 2 )dAa QCD {Dfizx, Q 2 )D^(z 2 , Q 2 ) + Df(z 2 , Q 2 )Df( Zl , Q 2 )} 



+ (contributions from Ad, As, Am, Ad and As) , (6) 
(iv) 

dAa^ n ~ 

~ Ag( V ,Q 2 )dAa PGF {^Df(z 1 ,Q 2 )Df(z 2 ,Q 2 ) + l -Df (z u Q 2 )Df (z 2 ,Q 2 ) 
+ l -Df{ Zl ,Q 2 )Df{z 2 ,Q 2 ) + - tt-(^)) } 

+^A M (r/, Q 2 )dAa QCD {D: { Zl , Q 2 )Df (z 2 , Q 2 ) + DJ" (* 2 , Q 2 )/^ {z x , Q 2 )} 

+ (contributions from Ad, As, Au, Ad and As) . (7) 

Due to the isospin symmetry and charge conjugation invariance of the fragmentation functions, 



various fragmentation functions in eqs.@-(0) can be classified into the following 4 functions, ||14j| 
D = D? = Df = Df = Df,D = Df = Df = Df = Df, D s = Df = Df = Df = 
Df and D g = Df = Df, where D and D are called favored and unfavored fragmentation 
functions, respectively. Considering the suppression of the s quark contribution to the pion 



production compared with the u and d quark contribution, we do not identify D s with -D.|15], [16 
By using these 4 kinds of pion fragmentation functions, we can make an interesting combination 
of cross sections which contains only the PGF contribution as follows; 

dAa n+ ^ + dAa n ~ n+ - dAa* +n+ - dAa*'*' (8) 
~ ^Ag( V ,Q 2 )dAa PGF {D( Zl ,Q 2 ) -D( Zl ,Q 2 )} {D(z 2 ,Q 2 ) -D(z 2 ,Q 2 )} . 

From this combination, we can calculate the double spin asymmetry All defined by 

dAa^ 71 ' + dAa^^ - dAa^^ - dAa^ 71 ' 



A 



LL 



da n w + da n n+ - da n n - da 7 " n 
Ag( V ,Q 2 ) dAapcF 

g(rj, Q 2 ) da PGF 

where the factor of the fragmentation function in eq. @ is dropped out from the numerator and 
the denominator of All-Q Therefore, from the measured All, one can get clear information of 
Ag/g with reliable calculation of dAapcF I dapcF- 



*For large Q 2 regions, heavy quarks might contribute to the PGF and QCD Compton. Even so, the final form 
of eq.(||) for All remains unchanged if Dq — Dq — Dq — Dq . 



3 Numerical calculation of cross section and double spin 
asymmetry 



Now, let us calculate numerically the cross section and double spin asymmetry An for the large- 
Pt pion pair production of pol-DIS. The spin-independent (spin-dependent) differential cross 
sections for producing hadrons hi and h 2 are given by|L7| 

''PGF 



d(A)a hlh2 d(A)a^ 2 - 



dz\d cos Q\d%2d cos 92dxdy dz\d cos B\dz2d cos 9 2 dxdy 



+ 3 , . (10 ) 



dz\d cos 9\dz2d cos 02dxdy 
Each term in the right hand side of eq. (JTTj|) is written as 

dz\d cos d\dz 2 d cos 8 2 dxdy 
= (A)g( V ,Q^C(6 1 ,92) 



dzid cos Q\dzjd cos 9 2 dx g dy 
x E _ eUD^(^Q 2 )D^Q 2 ) + (l~2)}, (11) 



d(A)£7. 



h\h,2 
QCD 



dz\d cos Q\dz2d cos d^dxdy 



q\ jjqyn^n j \ xj .£/ dZqd COS Q\dZgd cos 02dxqdy 



E _ e;(A)/ ff (i7,g 2 ) (7(^,08 



'QCD 



x{^(a,g 2 )^ 2 (6,Q 2 ) + (i - 2)} , (12) 

where 77 = x + (1 - x)tit 2 , Q 2 = zt/s, £1 = f 11 ^) zi, 6 = z 2 and C(0i,0 2 ) = 



^ — 2 ^ - — \ n ■ TTTn — TTTK 1 with Ti 2 — tan o- 

»7T1T2 8 cos 2 i 6»i cos 2 i6 2 sm^ i (61+6*2) ' ' 2 X ' Z 

Here we simply assume the scattering angle of outgoing hadrons to be the same with the 
one of scattered partons in the virtual photon-nucleon cm. frame. This assumption might not 
be unreasonable if observed particles are light hadrons with high energy, s is the total energy 
square of the lepton scattering off the nucleon. x, y and zi t 2 are familiar kinematic variables for 
semi-inclusive processes in DIS, defined as x = y — ^| and Zi >2 = P p^ 2 , where £, q, P 
and Pi 2 are the momentum of the incident lepton, virtual photon, target nucleon and outgoing 



hadrons, respectively. At the leading order of QCD, differential cross sections of hard scattering 
subprocesses, £g — > £ 
angle are given asjl8| 



subprocesses, £g £'qq and £ q — > £' 9 q , with two outgoing partons in an opposite azimuth 



d(A)a PGF (Q C D) 1 a 2 a s 7/(77 - x)(l - r]) 2 



dzid cos 9\d(j)\dziig)d cos 02dxg^ q )dyd(j)£ 1287T 2 (po ■ £)Q 2 x 

x B(e 1 ,e 2 )e 2 i e 2 \(A)M\ 2 PGF{QCD) , (13) 



with 

1 

sm(6 1 + 2 ) 

sin #1 -4- I Z7/ n (\ — T)\ A- (n — sin #o 

(14) 



x 



-77) + (77 - a;)}sin0 1 + {^(1 - 77) + (77 - x)}sin0 2 



sin 61 sin 9 2 

where Zi, z\ and z g are the momentum fraction of the outgoing parton i, i and g, respec- 
tively, to the incoming parton, and are given as Zi = ^ , zj^ = ^ . |I7| The amplitude 
|(A)M|p GF ^g CD ^ in eq.fllBD is presented elsewhere [|T5|] . 



By using these formulas and pion fragmentation functions ||15||, we have calculated the spin- 
dependent and spin-independent cross sections of the large-pr pion pair production and esti- 
mated the double spin asymmetry at the energy of COMPASS experiments.^] Here, we have 
taken the AAC|| and GS96[Q parameterizations at LO QCD as polarized parton distribution 
functions and GRV98JH and MRST98fl2lj as unpolarized ones. At = 13.7GeV, y = 0.75, 



Q 2 > lGeV 2 and W 2 > lOGeV 2 with kinematical values of 0i )2 and z 1>2 for the produced pion 
pair, the calculated results of the spin-independent (spin-dependent) differential cross sections 
and All are shown as a function of 77 in Figs.l and 2, respectively. In Fig.2, error was estimated 

by^ 8 All = + 3& where ^ = £t£ with ^ = + da(n~n+) and 

d(J2 = da(7i + 7i + ) + dcr(7T~7T~). L is the integrated luminosity and P = PsPrf, where Pb,Pt 
are the polarization of beam and target, respectively, and / is the dilution factor. Here we used 
L = 2 fb -1 (150 days running), Pb = 0.80 and Pp ■ f — 0.25. From Fig.2, one can see a big 
difference of the behavior of All depending on the models of Ag/g and hence, we can extract 
the behavior of Ag rather clearly from this analysis. 

4 Summary and discussion 

We proposed a new formula for extracting the polarized gluon distribution from the large-py 
light hadron pair production in pol-DIS by taking an appropriate combination of hadron pair 
productions. Since the double spin asymmetry All for this combination is directly proportional 
to Ag/g, the measurement of this quantity is quite promising for getting rather clear information 
on the polarized gluon distribution in the nucleon. 

The analysis can be applied also for the kaon pair or the proton pair production by considering 
the reflection symmetry along the V-spin axis, the isospin symmetry and charge conjugation 



t Calculation for HERMES energy was done|l8[ without error estimation. 




Figure 1: r] dependence of combined spin- independent and spin- dependent differential cross sec- 
tions defined at the denominator and numerator, respectively, of eq.(0) as a function of rj at 
y/s = 13.7 GeV, y = 0.75 for the deep inelastic regions (Q 2 > lGeV 2 and W 2 > 10 GeV 2 ) with 
kinematical values of 6^2, for the produced pion pair. 




Figure 2: 77 dependence of A LL at y/s = 13.7 GeV, y = 0.75 for kinematical values of lj2 , 2^1,2 for 
the produced pion pair. Solid line and dotted line are for AAC LO and GS96LO-C parameteri- 
zaton models, respectively. Ag/g itself is also presented for both parameterization models. 



invariance of the fragmentation functions. 

The same combinations of cross sections for light hadron pair production were discussed 
for photoproduction by other people |23[, while we studied the lepton-proton processes in deep 



inelastic region in this work. Both of those processes are useful for extracting Ag. 

This work is supported by the Grant-in-Aid for Science Research, Ministry of Education, 
Science and Culture, Japan (No. 11694081). 
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